Decreasing elliptic flow at the LHC: Calculations in a parton recombination approach 
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We apply the parton recombination approach to study the energy dependence of the elliptic flow, 
«2 in heavy ion collisions from AGS to LHC energies. The relevant input quantities (T, fi B ,Vt) at 
the various center of mass energies are obtained from fits to the available data. The model yields a 
good description of the integrated V2 data for charged particles at midrapidity from AGS to RHIC 
energies. In stark contrast to the current expectations, we observe a decrease of the integrated 
V2 values above the highest RHIC energy. Thus, we predict a decrease of vi at LHC energies 
compared to the RHIC results. This drop is attributed to negative «2 values for the underlying 
parton distributions at low to moderate transverse momenta that develops if the transverse flow 
velocity is high enough. At energies above the LHC regime, the present approach predicts even 
negative values for the integrated v%. 



The main goal of the current and past heavy ion pro- 
grams is the search for a new state of matter called the 
Quark-Gluon-Plasma (QGP) fjj. Major breakthroughs 
for the potential discovery [2J, |3| of this new state of mat- 
ter where (I) the experimental discovery of a large elliptic 
flow (i>2) at RHIC. This lead to the conclusion, that the 
matter formed in the early stage of a heavy ion reac- 
tion at RHIC behaves like a nearly perfect liquid (i.e. 
a liquid with a low viscosity) (0| and Refs. therein). 
(II) The observation of constituent quark number scal- 
ing w hadron( p hadron) = ^g^hadroB/^^ mea ning that 

the elliptic flow of baryons (three quarks, n q = 3) ver- 
sus mesons (two quarks, n q = 2)) scales like 3 : 2 at the 
hadron transverse momentum p^ adron . While the details 
of this scaling are still under discussion 0, @, 0, S, [§] , one 
might generally see this scaling as evidence for a recom- 
bination like hadronization process for the transition of 
partonic matter to hadronic matter []j| EH ■ 

In this letter, we will apply the recombination ap- 
proach to explore the energy dependence of the elliptic 
flow in massive (Pb+Pb/Au+Au) nucleus-nucleus reac- 
tions from the AGS energies to the highest LHC energy. 
For complemantory explorations of the elliptic flow at 
LHC within transport approaches the reader is referred 
to O, EE, Ll3l - Elliptic flow is a well chosen observable 
for the exploration of the energy dependence of flow ob- 
servables, because it exhibits a self quenching effect and 
is therefore mostly sensitive to the early (partonic) stage 
of the reaction, even at rather low beam energies. The 
structure of this letter is as follows: We start with a sum- 
mary of the energy dependence of the input parameters, 
then we discuss shortly the relevant recombination for- 
mulas, finally we present the results for the elliptic flow 
excitation function and the predictions for LHC. 

To apply this model to other energies than RHIC 
(y/s = 200 GeV) one has to model the dependence of 
the temperature T, the baryo-chemical potential /is and 
the transverse flow rapidity t\t — atanh(/?T) at the 
hadronization surface as a function of the center of mass 
energy. Details of the parameterisations can be found in 
the appendix. The values for RHIC energies and below 



are fitted to previously extracted flow velocities, at LHC 
energies, we obtain the values (for y / ijvw = 5.5 TeV): 
T = 175 MeV, fj,B ~= MeV and fh = 0.75c, which are 
in line with previous estimates [13]. 

With these parameters we predict the elliptic flow of 
hadrons using the collinear recombination approach fTTj |. 
Since we are presently only interested in bulk (i.e. low 
Pt quantities, we apply the recombination prescription 
in the low to intermediate pr-range from to about 
3 GeV at RHIC and about 5-6 GeV at LHC respectively, 
above this momentum, fragmentation will dominate the 
underlying quark distribution. Comparing the yields 
to experiment the formalism seems applicable down to 
Pt ~ 1 GeV at RHIC, which is approximately 3% of the 
matter produced, and we expect this lower limit to rise 
slightly to pr ~ 1.2 GeV at LHC. So using the recom- 
bination approach down to pt — > seems questionable, 
however, a part of the uncertainty in the recombination 
mechanism at low pt, introduced by the violation of en- 
ergy conservation, cancels after taking the ratios in Eq. 
([6]). Thus, the recombination formalism seems to give 
valid results for vo down to transverse momenta of several 
hundred MeV [11|. Further justification about the valid- 
ity of our results can be found in the paragraph about 
the mean V2- 

The picture that one quark and one anti-quark (or 
three quarks) with collinear momenta can recombine to 
form a meson (baryon) leads one to an integral over the 
product of the constituent-quark densities w a times the 
Wigner function of the hadron with the freeze-out 
hyper-surface £ and its normal vector u v fill ]: 
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d 3 N h 
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daP v u v {a) 

(27T) 5 
dx a 



(1) 



were the quark a carries the fraction x a of the hadron 
momentum P. 

We use the parameterisations from [ljj with a pure 
thermal spectrum for the recombining quarks and do not 
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FIG. 1: (Color online) Comparison of our asymmetry param- 
eter a to the eccentricity from CGC and the Glauber model 

El. 



include the power-law-tail which is only relevant for frag- 
mentation. At higher cm-energies recombination is ex- 
pected to be dominant up to even higher pr 03 , there- 
fore, we neglect the contributions from fragmentation for 
the present considerations. Following jTl[ we obtain the 
elliptic flow for non-central reactions from the asymmetry 
of the overlap region. To model peripheral collisions the 
asymmetry a of the collision depends on the geometric 
width and height of the transverse overlap zone 



h(b) = sjR\ - (6/2)2 w{b) =Ra _ {b/2) (2) 



as 



a 



h(b) - w (b) 
h(b) + w(b) 



(3) 



We have compared this simple asymmetry parameter 
with the eccentricity from other parameterisations in Fig. 

. One observes that the a parameter used (full line) is 
quantitatively close to the Glauber model (short dashed 
line) calculation for the relevant impact parameters. It 
has also been speculated that a much 'sharper' Color 
Glass Condensate initial distribution might be present in 
the initial state, this parametrisation is depicted by the 
long dashed line. The CGC values for the eccentricity are 
slightly higher and lead to a small quantitative change of 
the final V2 values, however the qualitative feature of a 
V2 decrease at high energies is not affected. 

This spatial asymmetry is translated into a velocity 
asymmetry in the parton phase by 



f».(^)=» J §,(l + a/( P r)coB(2^)) 



(4) 



this anisotropy less than slower ones we define 

1 



f(pr) = y 



(pt/po 



(5) 



with the parameter po — 1.1 GeV taken from [ll|. The 
elliptic flow is then given by 



v 2 = (cos(2$)) 



J cos(2$)f^d$ 



I 



d Pt 

d Pt 



(6) 



and is calculated for the individual quark species as fill ] 



v q 2 {p T ) 



f cos(20)/ 2 [a(0)] Ki [b{cf>)} d4> 



(7) 



with 



To account for the fact that faster partons will experience 



a(4>) =p T sinh(77 T (0))/T, 6(0) = m T cosher (<£))/T 

and the modified Bessel functions I n and K n . Replacing 
w a {R;p) in Eq.fQ with 

w' a (R;p)=w a (R;p){l + 2v 2 (p T )cas(2$)) (8) 

and again applying the definition for w 2 yields the 
hadronic elliptic flow. 

Let us start by comparing the transverse momentum 
dependence of the elliptic flow for pions, kaons, protons 
and omega-baryons in mid-peripheral Au+Au reactions 
at RHIC to the distributions obtained in mid-peripheral 
Pb+Pb reactions at the LHC as shown in Fig. [2j A 
first striking observation is that the V2(pt) values at each 
given pt are always lower at the LHC than at RHIC. A 
similar pattern was also observed within a parton trans- 
port approach, if the viscosity was set to the ADS/CFT 
limit [16j]. Even more surprising, we find strong out-of- 
plane emission of (multi-)strange particles, i.e. negative 
values of v 2 ■ These negative values of the elliptic flow pa- 
rameter for heavy particles have also been found in previ- 
ous exploratory studies and seem to be a general feature 
of the blast-wave like flow profile at high transverse veloc- 
ities [Isl . Iisl I20L l2l| . The negative V2 values have also be 
seen within a blast- wave fit of the proton elliptic flow to 
STAR data [121 and they have also been experimentally 
observed by the NA49 collaboration [23]. These results 
are supported by the recent preliminary data on J/r/) el- 
liptic flow from PHENIX [lB|. It reflects the depletion 
of the low pr particle abundance, when the source ele- 
ments are highly boosted in the transverse direction. At 
higher pr and or lower transverse flow velocities, the op- 
posite effect dominates and results in a more pronounced 
emission in-plane again. The details depend on the un- 
derlying transverse flow parameterisation as discussed in 
[2l| . One might argue that negative elliptic flow values 
are an artefact of the blast wave parametrisation, how- 
ever, also transport simulations indicate slightly negative 
v 2 values for heavy particles at very low pr [23] ■ Thus, 
the qualitative behaviour of a negative V2 is a well-known 
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FIG. 2: (Color online) Elliptic flow parameter «2 as a function 
of the transverse momentum pt for mid-central (6 = 8 fm) 
Au+Au/Pb+Pb reactions as obtained from the recombina- 
tion approach. The lines indicate the calculations for for 
7r + , K + ,p and fi~. Upper Panel: Vi at RHIC (y/s = 
200 GeV) compared to PHENIX data Q3]. Lower Panel: 
Prediction for LHC (y/s = 5.5 TeV). 



observation. The surprising observation in the present 
work is the quantitative strength of this effect at LHC. 

From these arguments we conclude that the blast-wave 
flow profile is responsible for the negative V2 on the quark 
level, which then enters the hadron elliptic flow via the 
parton recombination. To understand this effect in more 
detail we go to the thermal quark-spectrum [ll[ with the 



energy 

E(R,p) — rriT cosh(r] — y) cosh t]t ((f)) 

—pr cos(4> — $) sinh?7<r(0) (9) 

For simplicity we look at midrapidity (y = 0) and for 
high r\ the spectrum is very low, so we consider the re- 
gion around r\ — 0. For high cm-energies, when the 
source is highly boosted transversally, the particles will 
mainly be emitted in the direction in which the fireball 
flies, so we can simplify even more and set (j> = Be- 
cause we have no longitudinal momentum we can re- 
place the momentum with the transverse rapidity of 
the parton: pr = msinhr^ and hit = y rn? + pip = 

myl + sinh 2 rj T — m cosh T]f . So the energy of the quark 
in the transverse moving source is 

E = mcosh [rj T (l + af(p T ) cos(20)) - r) T ] (10) 

For a fixed quark rapidity r\ T < rj T (low pr) the energy of 
the quarks emitted in plane is higher than the energy of 
quarks emitted out of plane. With the thermal spectrum 
more energy means less particles, therefore a negative 
t>2. At r\ T = i] T one would expect the zero-crossing, and 
above a positive V2. 

In this way one can visualise our predicted negative 
elliptic flow for low pr at LHC. At RHIC this picture 
can not be applied, because the simplification = $ is 
only valid for high 77^ (high cm-energies). 

A different way to obtain an analytic expression which 
explains the negative V2 is to consider only the in-plane 
((f) = 0) and out-of-plane (6 = ir/2) directions (similar to 
the analysis performed in [l9(]). The <f) integration breaks 
down and the elliptic flow is then given by 

h {a(0)] K x [6(0)] - h Ntt/2)] K x [b(n/2)} 
V2 {PT> I [a(0)} K x [6(0)] + I [o(tt/2)] K x [b(w/2)] 

(11) 

For pr — » the argument of the Bessel functions I n goes 
to zero and /„ becomes constant with I2 — * and Iq — > 1, 
Therefore they are independent of the angle. This leads 
to 





lim PT ^ vl(p T ) 


_ -Ma(O)] 
Ja[o(0)] 




Ki 


m cosh(?7?,(l+a)))/T] -K i 


m cosh(r;^ (1 


-o)))/T 


Ki 


m cosh(?7° (l+a)))/Tj+A'i 


rn cosh(r/S^ (1 


-a)))/T 



(12) 



Since K i is a monotonically decreasing function the nom- 
inator the elliptic flow is negative (for some small trans- 
verse momenta). The specific values depend on the mass, 
the mean flow rapidity 7j T and the temperature. For 
increasing mass, or increasing transverse flow rapidity 
or decreasing temperature, the elliptic flow will become 
more negative. 

Folding the distributions for t>2 (pt) with the corre- 
sponding transverse momentum distributions yields the 
integrated v 2 at midrapidity for each particle species, see 
Fig. |(3l). A first observation is the apparent good de- 
scription of the available charged particle v 2 data from 
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FIG. 3: (Color online) Comparison of {V2) for charged hadrons 
at midrapidity as a function of center of mass energy from 
the present recombination approach for Au+Au/Pb+Pb re- 
actions at an impact parameter of b = 5 — 6 fm. Data points 
and the extrapolation to LHC are taken from [2(|. 



the AGS to the RHIC energy regime. However, when 
going above the highest RHIC energy, we predict that 
the integrated elliptic flow saturates and then starts to 
decrease in the LHC energy range. This finding is in 
stark contrast to the current expectations in the heavy 
ion community - compare e.g. to the linear extrapola- 
tions towards LHC by Borghini and Wiedemann [26]). 

A rather critical assumption in this study is the ap- 
plicability of the recombination approach for the ellip- 
tic flow for small transverse momenta on the order of 
pr < 1 GeV. We want to emphasise that our re- 
sult of the decreasing mean elliptic flow (va) at LHC is 
not affected by the validity of this assumption, because 
(v 2 ) ~ V2 {(pr)) and (pr) > 1 GeV in LHC regime. How- 
ever, to show the robustness of the prediction Fig. (@]) 
depicts the elliptic flow at a fixed pr as a function of yfs. 
With pt = 1, 1.5 and 2 GeV the elliptic flow exhibits the 
same drop as the mean u 2 from Fig. J3]). 

We have calculated the energy dependence and trans- 
verse momentum dependence of the elliptic flow param- 
eter V2 for mid-central Au+Au/Pb+Pb collisions from 
AGS to LHC energies within a parton recombination ap- 
proach. We find a reasonable description of the V2 data 
over the whole inspected energy regime, indicating that 
the measured V2 values are consistent with the assump- 
tion that a major part of the elliptic flow was created 
in the partonic stage. We predict that the integrated 
V2 of charged particles at midrapidity will decrease from 
RHIC to LHC energies, due to the strong transverse flow. 
In detail, we link this to the prediction of a negative v 2 
component developing at low transverse momenta to the 
blast-wave like flow profile of the underlying quark distri- 



p =1 GeV 

p T =1.5GeV - 



p T =2 GeV - 



H 
Q. 

"O 

CD 

X 



crj 

(M 
> 




-cm 



FIG. 4: (Color online) Elliptic flow i>2 of charged hadrons 
at b = 5 — 6 and fixed pr as a function of -^/sjvjv for pr = 
1, 1.5 and 2 GeV. 



bution. Because this effect is strongest for heavy quarks 
it most visible for (multi-) strange hadrons. Above the 
presently envisaged LHC energy for nucleus-nucleus re- 
actions, we predict that the mean elliptic flow will even 
turn negative. It should be pointed out that the present 
prediction is in striking contrast to all former assump- 
tions about the behaviour of t>2 at LHC. 
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APPENDIX A 

To extract T and /is at the phase boundary we take 
the freeze-out values of the baryo-chemical potential from 
different accelerators and fit them against their cm- 
energy via 



Freeze 

Ms — 



1 + V S NN /c 



(13) 



with a = 1.478 GeV, b = 0.802, c = 2.096 GeV. Then 
we follow the isentropes from freeze-out to the phase 
boundary in order to connect this /i^ reozc to a \ib value 
at the phase boundary (Fig. J5])). In the relevant region, 
the connection is approximately linear with 



Hb = d 



Freeze 
PB ) 



d = 0.938 



(14) 



Due to the lack of lattice data at high baryo-chemical 
potentials, we calculate the phase transition line with 
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PIG. 5: (Color online) The QCD phase diagram with a freeze- 
out curve fitting the collider experiments (RHIC, SPS, SIS, 
AGS) and some isentropes to the MIT-Bad model phase 
boundary. 



a simple MIT-Bag model with a critical temperature of 
Tc = 175 MeV at vanishing chemical potential. The line 
of the phase transitions temperature is 



T 



\ 




with C 



/M£\ 4 9q_ 
\ 3 J C 2 



9q 



c_ 

72 



(15) 

16. Details like 



^(7g q + 4g g ), g q = 12, g g 
the order and nature of the parton-hadron transition are 
not relevant for the present study. Because the fugacity 
7 = exp(/xs/T) does not enter in the equation for V2 (in 
any case 7 ~ 1 for all relevant flavours at LHC energies), 
we neglect it in this paper. The temperature, baryo- 
chemical potential and the fugacity at hadronization are 
shown in Fig. (J6j) . 



APPENDIX B 

While the hadronization temperature shows only a 
weak dependence on energy, the transverse expansion ra- 
pidity t\t is strongly energy dependent. We fit the flow 
rapidities extracted from the experimental data [13] at 
kinetic freeze-out. The data given in [28[ is for mean 
transverse flow velocity (3t, but we use the rapidity 
7]t = tanh j3r to assure that the velocity stays less than 
c. To fit these values we choose 

^ eeze (Vi) = a + bx + cx 2 + d\n(x) ,x = ln (y^) (16) 

with the constants a = 0.418, b = —0.064, c = 
0.012, d = 0.170. As these values are extracted at freeze- 



out, we scale the obtained transverse rapidities by a con- 
stant factor k = 0.85 to obtain the transverse flow at 
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FIG. 6: (Color online) Temperature T and baryo-chemical 
potential fiB (left y-axis) and fugacity 7 = e ' MB / T ' (right y- 
axis) at the phase transition as a function of the center of 
mass energy ^/snn. 
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FIG. 7: (Color online) Parameterisation of tjt (full line) as 
function of a/s . The data at kinetic freeze-out (crosses) are 
taken from 1 28]. 



the hadronization surface. Using these parameters our 
value for vt = 0.54 at RHIC energies agrees with the 
value from [llj. For the LHC energy (s/s — 5.5 TeV) 
we obtain a transverse flow velocity of vt = 0.75 also in 
line with previous estimates [13]. Fig. J7]) depicts our fit 
(line) and the available data on rj T (crosses). 
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